First-order transition between a small-gap semiconductor and a ferromagnetic metal 

in the isoelectronic alloys FeSii-xGe^. 
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The contrasting groundstates of isoelectronic and isostructural FeSi and FeGe can be explained 
within an extended local density approximation scheme (LDA+U) by an appropriate choice of the 
onsite Coulomb repulsion, U on the Fe-sites. A minimal two-band model with interband interactions 
allows us to obtain a phase diagram for the alloys FeSii-^Gea;. Treating the model in a mean field 
approximation, gives a first order transition between a small-gap semiconductor and a ferromagnetic 
metal as a function of magnetic field, temperature, and concentration, x. Unusually the transition 
from metal to insulator is driven by broadening, not narrowing, the bands and it is the metallic 
state that shows magnetic order. 



The unusual magnetic susceptibility of FeSi has been 
a subject of interest for a long time. Jaccarino et 
al.u found a rapid crossover around room temperature 
from activated behavior at low temperature to an ap- 
parent localized Curie- Weiss form at high temperature. 
This has stimulated a number of proposed explanations. 
Among them are a form of Kondo insulating behavior in 
FeSrra, although an underlying microscopic model has 
not been clarified, and almost ferromagnetic semiconduc- 
tor modelcra. The latter received support from the de- 
tailed LDA band structure calculation of Mattheiss and 
Hamanncl. They found a small gap semiconductor with 
formal neutral valences Fe 0+ , Si 0+ on both elements. 
Later Anisimov et all! extended these calculations to 
include the onsite Coulomb interaction, U, in a.-mean 
field approximation, the so-called LDA+U methodEl, and 
found, for a reasonable choice of U, a ferromagnetic 
metallic phase very close by in energy. This led them 
to propose that a transition to this ferromagnetic state 
could be driven by applying a magnetic field (B) and 
that the unusual temperature (T) dependent suscepti- 
bility in FeSi reflected the proximity to a critical point 
of this transition at finite (B c , T c ). However the pre- 
dicted value of B c («170T) is too large to be reached 
in the laboratory and the existence of the critical point 
remains untested experimentally. A non-trivial predic- 
tion was the fractional value of the saturation moment of 
the ferromagnetic phase at S = 1/2 or 1 /is/Fe. This is 
a consequence of the band structure and cannot be rec- 
onciled with a local model and a 3d 8 configuration for 
Fe°+. 



the fractional S — 1/2 value quoted above EJO. In this 
letter we report nLDA+U calculations (in TBLMTO cal- 
culation schemco) for FeGe and the isoelectronic alloys 
FeSii_a;Ge x . Motivated by these results we use a sim- 
plified phenomcnological model to explore the complete 
phase diagram of the isoelectronic alloys FeSii_ x Ge x 
with varying T and B. Alloying FeSi with Ge allows 
one to tune the predicted semiconductor to ferromag- 
netic metal transition and its unusual critical point at 
(B c , T c ), to experimentally convenient values. 
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FIG. 1: The density of states (DOS) near the Fermi level at 
zero energy obtained from LDA calculations. The solid and 
dashed lines represent FeSi and FeGe with energy gaps 0.08 
eV and 0.03 eV respectively. The corresponding widths of the 
peaks above the Fermi level are 0.5 eV and 0.37 eV . 



The magnetic properties of the isoelectronic and 
isostructural FeGe compound have aroused less inter- 
est. It is a magnetic metal with a long period spi- 
ral form which is simply a ferromagnet twisted by the 
Dzyaloshinskii-Moriya interaction that is a consequence 
of the absence otinversion symmetry at the Fe site in this 
cubic structure£3. Interestingly the saturation moment is 



Both FeSi and FeGe crystallize in the cubic B20 struc- 
ture which can be viewed as highly distorted rocksalt 
with 4 FeSi formula units in the primitive cell. Each 
Fe site has 7 Si neighbors and point group symmetry, 
C3. The LDA band structure calculations of Mattheiss 
and Hamannu place the Fermi energy within a narrow 
manifold of 20 bands of predominantly Fe 3d charac- 
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ter lying within a larger hybridization gap with the Si 
(3p, 3s) states. A fuller description of the band structure 
will be presented elsewheret2l. A nontrivial feature of 
their results is the presence of a small but complete band 
gap separating 16 filled valence bands from 4 empty con- 
duction bands. The conduction bands have nonbonding 
character with respect to the Si atoms and are divided 
by a pseudogap into two lower lying quite narrow bands 
and two higher lying bands which are much wider. The 
results of our LDA band calculations for FeGe are similar 
but with an even smaller band gap and narrower overall 
bandwidths as expected from the larger lattice constant. 
The total density of states for both compounds is shown 
in Fig. 1. 

Anisimov et alE extended the LDA calculations for 
FeSi to incorporate the onsite Coulomb repulsion (U) on 
the Fe-siteS|4n a mean field approximation scheme known 
as LDA+UD. They found a local minimum in the total 
energy versus magnetization at a value « 1/x^/Fe which 
moved to lower energy with increasing U. In this magne- 
tized state the lower pair of conduction bands are filled 
for the majority spin direction while the Fermi level of 
the minority spin electrons lies in the valence band com- 
plex leading to metallic behavior. The relative positions 
of the two minima corresponding to the small gap semi- 
conductor and magnetized metal, depends on the choice 
of U which is not known precisely. 
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FIG. 2: The evolution of the total energy ( with energy of 
nonmagnetic solution taken a zero ) as a function of the spin 
moment M(fj,B / Fe) for the value of U = 3.7eV. The solid, 
dashed and dashed-dotted lines correspond to FeGe, FeSi and 
FeSio.58Geo.42 respectively. 

We have used the TBLMTO scheme! to perform 
LDA+U calculations for FeGe. The narrower band- 
widths and energy gap lead to a lower value for the crit- 
ical U, where the two energy minima cross, in FeGe rel- 
ative to FeSi. Thus it is possible to obtain within the 
LDA+U scheme for a common sialue of U, the correct 
groundstates for FeSi and FeGeEil (see Fig. 2). There 
are two unusual aspects to this metal-insulator transition 
(MIT). It is the metallic state, rather than the insulating 



state, which is magnetically ordered. Secondly, the tran- 
sition from insulating to metallic behavior is driven by 
narrowing, rather than increasing, the bandwidth. The 
inverted nature of the metal-insulator transition is a di- 
rect consequence of the fact that the MIT is driven by 
the paramagnet-ferromagnet transition. 

FeSi and FeGe are the end members of the isoelectronic 
and isostructural alloys, FeSii-zGez-i-We have extended 
the LDA+U calculations to the alloysEj using experimen- 
tal lattice parameterala. The critical value x c of the first 
order transition is sensitive to the choice of U. For the 
value of U = 3.7eV, illustrated in Fig. 2, x c = 0.4 in_A 
good agreement with the experimental value x c — 0.3lj. 

To proceed further we introduce and solve (within a 
mean field approximation) a minimal phenomenological 
model for the isoelectronic alloys FeSii-^Gej;. Motivated 
by the experimental fact thaLpthe ordered moment in 
FeGe is » per Fe atornOtJ which implies a frac- 
tional magnetization with respect to the paramagnetic 
Fe 3(i 8 configuration, we are led to an itinerant model 
of magnetism. We consider the following model for the 
conduction and valence bands. The valence band, de- 
scribing the upper part of the manifold of occupied Fe 
3d bands, stretching below zero energy in Fig. 1, has a 
width 2W and contains one state per spin per Fe atom. 
The conduction band (of width W) contains 1/2-state 
per spin per Fe atom and when multiplied by 4 (which 
is the number of Fe atoms per primitive cell), it corre- 
sponds to the two narrow conduction bands just above 
the Fermi level. Both the valence and conduction bands 
are assumed for simplicity to have a constant density of 
states N(e) — (2W)~ 1 per spin, with a k-independent 
gap of 2A. Taking 2 conduction electrons per Fe atom, 
leads to a fully occupied valence band and a semicon- 
ducting ground state in the noninteracting case. 
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FIG. 3: Density of states of the model band structure, (a) 
Nonmagnetic semiconducting state, (b) Magnetic metallic 
state. The arrows note majority and minority spins. 

We propose that three types of effective local interac- 
tions should control the physics: intraband repulsion U, 
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interband repulsion V, and exchange coupling between 
the bands /. The interactions U, V, and I are combina- 
tions of the matrix elements of the Coulomb interaction 
between the relevant Wannier orbitals and are taken as 
a phenomenological input to our theory. Furthermore, 
we assume that the primary effect of the U term is to 
strongly renormalize the bandwidth W with respect to 
estimates from the LDA calculations. In order to de- 
scribe the experimental data, we have therefore decided 
to use a phenomenological bandwidth rather than the 
calculated W ~ 0.5 eVQ. Thus, in order not to double 
count its effects, the U term is not included in our model 
Hamiltonian which reads 



H = 



A,k,c 



(e Ak - <JH B B)n X ka - -= ^ (lm R + Vp^) 



R 



where A = — 1 , + 1 denotes the valence and conduction 
bands, respectively, and £Ak and nAko- is the band energy 
and number of electrons with momentum k and spin a 
in band A. The operators /5r, = 2 + J2\a ^ n ARcr and 
ifiR = X^Aer ^"ARtT measure the local number of charge 
carriers and the local magnetization, respectively. The 
spin quantization axis has been chosen in the direction 
of the external magnetic field B. The Hamiltonian can be 
written also in an explicitly spin-rotation invariant form, 
which we do not specify here. 

Below we describe the phase diagram of the 
FeSii-^Gea; alloys by fixing the values of /, V, and A 
irrespective of the Ge content x. Alloying changes the 
bandwidth W which decreases with increasing x, in ac- 
cordance with the LDA results. Treating the interaction 
terms in the mean field approximation allows to replace 
the operators pn and ttir, by the expectation values p 
and m, respectively. The main technical difference of the 
present study with respect to Ref. ^ is that our model is 
not particle-hole symmetric, and therefore in addition to 
m and p, the chemical potential p has to be calculated 
sclf-consistently. Assuming finite values of m and p, the 
single-particle energies in the lower and upper bands are 
LO\, a = e — XVp — <r(pBB + Im). The corresponding mean 
particle numbers are 



(n_, CT ) = (2Wy 1 
(n +<<T ) = (2Wy 1 



-A 



de 



-A-2W 1 + exp(w_ CT - p)/T 
A+w de 

A l + exp(u> +t(T — pL)/T 



, (1) 

(2) 



The equation for the total number of particles 2 = 
J2\ a (n\,a-) together with Eqs. form a closed set 

of equations for p, m, and p, which have to be solved in 
general numerically. 

At T = 0, the model can be solved analytically. For 
sufficiently large bandwidth W, the following two phases 
compete: (i) semiconductor with p = 0, m — and (ii) 
ferromagnetic metal with p = 1, m = 1. Their ener- 
gies are E semi = -2A - 2W and E FM = -A - 3W/2 - 



(I + V)/2 — pbB, respectively. With changing band- 
width, at B = there is a level crossing between the 
semiconducting phase and the ferromagnetic metal at 
W = W c (0) = I + V - 2A. For W > W c (0), the fer- 
romagnetic phase is stable only in magnetic fields larger 
than a critical magnetic field pbB c (0) — (W — W c (0))/2. 
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FIG. 4: Phase diagram in the W-T plane at B = 0. The 
semiconductor-ferromagnet transition (solid line) is of first or- 
der at low temperatures, up to a critical point at W c ~ 104 
meV and T c ~ 37 meV. At higher temperatures, the transi- 
tion is of second order (dashed line) . The spin susceptibility in 
the paramagnetic phase peaks at temperatures indicated by 
the dash-dotted line. The inset shows the magnetization as a 
function of temperature for a nearly critical alloy composition 
in an applied magnetic field. 

We take A = 20 meV in qualitative agreement with ex- 
periment and the effective interaction parameters / = 80 
meV and V = 65 meV close to the values proposed in 
Ref. 0. The resulting phase diagram for B — shown in 
Fig. 4 implies that the phenomenological bandwidth of 
FeGe, W FcGc is less than W c {0) (W c (0)=105 meV). Note 
that the Curie temperatures predicted by our model for 
W < W c (0), T c ~ 400 K, are in reasonable agreement 
with the experimentally observed transition .temperature 
to a long-period spiral in FeGe, T N w 280 KM. The phe- 
nomenological bandwidth of FeSi, WFcSi7 w hich is greater 
than W c (0), can be determined by fitting the spin suscep- 
tibility. For W > W c (0), the spin susceptibility develops 
a strong peak around room temperature, confirming the 
interpretation of Anisimov et al. The peak position as a 
function of W is shown in Fig. 4. Fitting the peak value 
of x to w 27p 2 B /e\W, we estimate W FcS i « 130 meV. 

Nearly critical semiconducting alloys with W slightly 
in excess of Wc(0) should exhibit metamagnetic transi- 
tions to metallic phases at experimentally accessible mag- 
netic fields. For instance in a sample with W = 110 meV, 
the transition occurs at B c (0) = 43.1 T at T = 0. The 
temperature evolution of the magnetization m for such a 
sample for magnetic fields in the vicinity of the first-order 
transition is shown in the inset in Fig. 4. Note that at 
low temperatures the critical magnetic field lowers with 
increasing temperature. This follows simply from the 
larger electronic entropy of the metallic phase. At higher 
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temperatures the phase boundary in our model bends 
towards smaller bandwidths W , in qualitative agreement 
with the higher susceptibility of small bandwidth systems 
to various forms of symmetry breaking transitions. The 
different low- and high-T slopes of the phase boundary 
lead to a reentrant B = phase diagram and also to reen- 
trant metamagnetic transitions in a narrow field range, 
as shown explicitly in the inset in Fig. 4 for B = 25.9 T. 

We should like to point out that not only the magnetic 
properties are anomalous at the metamagnetic transition. 
Since the transition is between a semiconductor and a 
ferromagnetic metal, large magnetoresistance is to be ex- 
pected. As a result the critical endpoint of this meta- 
magnetic transition should show particularily interest- 
ing magnetoresistance behavior at temperatures around 
room temperature. 



diet that, for intermediate values of the Ge content x, this 
transition can be realized as a metamagnetic transition 
at experimentally accessible magnetic fields. Large mag- 
netoresistance is predicted at such a transition. Another 
possible way to realize the transition at experimentally 
accessible magnetic fields would be to apply pressure to 
FeGe. 
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FIG. 5: Complete phase diagram of the isoelectronic al- 
loys FeSii-zGe^. The semiconducting phase at large band- 
widths W (small Ge content x) is separated from the ferro- 
magnetic metal phase at small W (large x) by a surface of 
first-order transitions terminating in a line of critical points 
(solid squares) around the room temperature. At B = and 
at Curie temperatures which lie above the temperature of the 
critical point, there is a second-order phase transition between 
the paramagnetic and ferromagnetic phases (open circles). 

In conclusion, we have calculated the complete phase 
diagram of the isoelectronic alloys FeSii^^Gej, with vary- 
ing T and B. Our main results are summarized in a 
three-dimensional phase diagram, Fig. 5, which shows 
that, depending on x, T, and B, the alloy FeSii-^Ge^ 
can be a semiconductor or a ferromagnetic metal. The 
two phases are separated by a surface of first-order tran- 
sitions which terminates at high temperatures in a line of 
critical points. Note, this_phase diagram resembles that 
found by Pfleiderer et aO for isostructural MnSi under 
pressure and magnetic field but with the important dif- 
ference that all phases of MnSi are metallic so that a 
semiconductor-metal transistion is not involved. We pre- 



1 V. Jaccarino et al., Phys. Rev. 160, 476 (1967). 

2 T. Mason et al, Phys. Rev. Lett. 69, 490 (1992). 

3 D. Mandrus et al, Phys. Rev. B 51, 4763 (1995). 

4 Y. Takahashi and T. Moriya, J. Phys. Soc. Jpn. 46, 1451 
(1979). 

5 S. N. Evangelou and D. M. Edwards, J. Phys. C 16, 2121 

(1983) . 

6 L. F. Mattheiss and D. R. Hamann, Phys. Rev. B 47, 13 
114 (1993). 

7 V. I. Anisimov, S. Yu. Ezhov, I. S. Elfimov, I. V. Solovyev, 
and T. M. Rice, Phys. Rev. Lett. 76, 1735 (1996). 

8 V.I. Anisimov, J. Zaanen and O.K. Andersen, Phys. Rev. 
B 44, 943 (1991); Vladimir I. Anisimov, F. Aryasetiawan, 
A.I. Lichtenstein, J. Phys.: Condens. Matter 9, 767 (1997). 

9 O.K. Andersen, and O. Jepsen, Phys. Rev. Lett 53, 2571 

(1984) . O.K. Andersen, Z. Pawlowska, and O. Jepsen, 
Phys. Rev. B 34, 5253 (1986) 

B. Lebech, J. Bernhard, and T. Freltoft, J. Phys.: Condens. 
Matter 1, 6105 (1989). 

1 R. Wappling and L. Haggstrom, Physics Lett. 28A, 173 
(1968). 

2 L. Lundgren, K. A. Blom, and O. Beckman, Physics Lett. 
28A, 175 (1968). 

3 V.I. Anisimov et al (to be published). 

4 Because of the neglect of spin-orbit coupling in these calcu- 
lations, the Dzyaloshinskii-Moria interactions which cause 
the long period spiral order in FeGe, are not included and 
ferromagnetic order results. 

5 The electronic structure of FeSii-^Ge^ alloys was calcu- 
lated in-jthe approximation, where in TBLMTO calculation 
schemetl the potential parameters for non-metallic atoms 
were taken as the corresponding average of the Si and Ge 
parameters. 

6 Z.Fisk (private communication) 

7 C. Pfleiderer, S.R.Julian and G.G.Lonzarich, Nature 414, 
427 (2001) 



4 



